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Intense (~ 20 dB) enhancements of the incoherent radar backscatter spectrum from 
the topside ionosphere have been observed with the Millstone Hill UHF radar. En- 
hancements occurring at the local ion acoustic frequency causing large asymmetries 
in the measured ion line may be produced by current-driven instabilities ( Foster et a/., 
1988] (henceforth, FOS88). These enhancements pose a practical problem for space 
surveillance systems because their cross section and spectral width are characteristic 
of satellites. 

Conversely, their hard target signature complicates the study of naturally occur- 
ring ERB events; it is nearly impossible to distinguish them from satellites based on a 
single measurement. Statistical comparisons of observed coherent echo distributions 
with predictions from a satellite catalog have been used to broadly identify periods of 
ERB activity. A series of experiments using multiple diagnostics, including satellite 
instruments, for simultaneous observations have established the association of ERB 
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with large fluxes of soft suprathermal electrons carrying field-aligned currents. Zenith 
data are also presented which show the asymmetric growth of ion acoustic waves di- 
rectly above Millstone Hill. Zenith measurements filter out satellite contamination 
because the geometry constrains the Doppler shift from orbiting bodies to velocities 
significantly less than the ion sound speed in the F region. 

Details of these results are presented in this chapter. First, however, a brief 
description of the normal incoherent scatter spectrum is reviewed, followed by a de- 
scription of the characteristics of the ERB spectra. 


5.1 Incoherent (Thomson) Scatter Process 

The radar technique for studying the earth’s upper atmosphere typically employs a 
(HF, VHP or UHF) radar to transmit a powerful pulse of electromagnetic radiation 
into the atmosphere. Following transmission, the radar is used to receive the radiation 
scattered by the atmosphere back to the antenna. If the details of the scattering 
process are understood, meaningful information about the medium can be extracted 
through the appropriate processing of the weak received signal. For the partially 
ionized upper atmosphere above about 60 km altitude, the electromagnetic waves 
are scattered by electrons in a process known as incoherent or Thomson scatter, 
named after J. J. Thomson who first showed that individual electrons can scatter 
electromagnetic radiation [1906]. 

5.1.1 Relationship Between Density Fluctuations and Scat- 
tering Properties 

The frequency spectrum of the high frequency radiation scattered by the ionosphere 
is directly related to the power spectrum of the electron density fluctuations in the 
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plasma. Mathematically the cross section of the plasma as a function of radar carrier 
offset frequency, a;, may be expressed as 

= <r.V(|»(k,«)| a ) (5.1) 

Here, <r e ~ 10 -28 m 2 is the scattering cross section of a single electron, V is the 
scattering volume, ( ) indicate a time average; k and u> are the difference of the 
incident and scattered wave vectors and frequencies, respectively, viz. k = k,- — k, 
and u = u>i— u> 3 . For the usual backscatter radar geometry, k = 2k*; thus, the incident 
radar wavelength, 27r/|k*|, determines the scale size of density fluctuations sampled. 
Here, n(k,w) is the space-time Fourier transform of the usual number density in the 
plasma, n(r, £), defined by 

1 r ^t/2 

n(k,u) = T7 = dr dt n(r, t) exp[ik • r - iu>t] (5.2) 

V 1 Jv J-T/2 

T in this equation represents a time greater than the correlation time of the plasma 
medium. The total cross section of the plasma can be obtained by integrating (5.1) 
over uj, 

°tot = / ^ dw (5.3) 

J — oo 

Spectrum of Density Fluctuations 

The key to determining the characteristics of the received scattered wave, then, lies 
in calculating the power spectrum of density fluctuations in the space plasma. 

When only the random thermal motions of the scattering particles are considered, 
the spectrum exhibits a broad peak centered on the carrier frequency and a character- 
istic width of the electron thermal velocity, A uj ^ kv te [ Fejer , I960]. The total cross 
section as determined by (5.3) is found to be simply N<r e , where N represents the 
total number of electrons in the scattering volume. This result has a straightforward 
physical interpretation as follows. 
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The electrons are assumed to undergo random thermal motions, implying that 
each electron will scatter signals of random phase relative to the other electrons. At 
the radar receiver the signal powers rather than the signal voltages will add, and the 
resulting cross section will be equal to the sum of the individual contributions of each 
electron, N<r e . The random phase of the scattered signals suggests the use of the 
name incoherent scatter to describe this process. 

However, as more detailed calculations show, the presence of ions in the plasma 
does introduce a degree of coherence between the electrons’ motions [e.g., Hagfors, 
1961; Salpeter , 1960, 1961; Dougherty and Farley , 1960; and others]. The influence of 
the ions is found to be important when the incident radar wavelength is much larger 
than the so-called Debye length in the plasma, given by 



where e a is the permittivity of free space, K is Boltzmann’s constant, and T e , n, e, 
v te , and u p are the electron temperature, density, charge, thermal velocity and plasma 
frequency, respectively. 

The Debye length is measure of the distance over which the plasma can shield 
out electric fields; in the ionospheric F region, ze .003 meters. For distances 
greater than this length the plasma dynamics are characterized by collective processes 
rather than the random thermal interactions of individual electrons. The resulting 
scattering can best be thought of as arising from density fluctuations associated with 
electrostatic oscillations in the plasma 1 . The idealized power spectrum of density 
fluctuations under these conditions is shown in Figure 5-1. 


l The term incoherent is not strictly accurate in this case and an alternative name, Thomson 
scatter , has also been applied to this phenomenon. 
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Figure 5-1: Incoherent scatter power spectrum when A » A d- [after Rishbeth, 1989] 

5.1.2 The Incoherent Scatter Radar Spectrum 

The power spectrum in Figure 5-1 consists of two components. The double— humped 
centred spectrum results from density fluctuations associated with ion waves and 
is correspondingly known as the ion line. The sharp peaks flanking the ion line 
are caused by high frequency electron density fluctuations called plasma waves; this 
component is therefore known as the plasma line. 


The Ion Line 


The ion waves that give rise to the ion line in the incoherent scatter spectrum are 
longitudinal oscillations with a phase velocity given by the “sound” speed in the 
plasma. The dispersion relation for these ion acoustic waves is given here, 


w (KT e + KT\ \ l/2 
k \ M; ) 


(5.5) 


where Mi is the ion mass and C s is defined as the ion sound speed in the plasma, typ- 
ically 1-2 km/sec in the ionospheric F region. While these waves are ion oscillations, 
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the electrons follow the motion because of their electrostatic attraction to the massive 
ions. The two peaks in the ion line power spectrum may be thought of as reflections 
from ion acoustic waves travelling towards (upshifted peak) and away (downshifted 
peak) from the radar. The peaks are broadened due to ion thermal motions. Analy- 
sis of the shape and magnitude of the spectrum yields several important ionospheric 
parameters [see, e.g., Evans , 1969]. 

The spectral width corresponds to the ion acoustic phase velocity, which from 
Equation (5.5) depends primarily on the electron temperature and the mass of the 
dominant ion species. The offset, Sfi, of the entire spectrum relative to the radar 
center frequency yields the line-of-sight bulk plasma drift. The electron- ion temper- 
ature ratio determines the sharpness of the spectrum’s peaks. The area under the ion 
line curve is directly related to the electron density. These parameters are usually 
derived through a numerical model which adjusts their values to achieve a best fit to 
the observed power spectrum. 

The Plasma Line 

The plasma line component of the incoherent scatter spectrum is usually much weaker 
than the ion line and more difficult to measure. The sharp peaks result from electron 
oscillations near the plasma frequency. The dispersion relation for these waves in a 
magnetized plasma is given by, 

+ + (5.6) 

where Cl e is the electron cyclotron frequency and 9 is the angle between the wave 
and the geomagnetic field. Usually the dominant term in (5.6), the electron plasma 
frequency, uj p oc n 1 / 2 , provides a very accurate indicator of the electron density. The 
intensity of the plasma line is dependent on the electron temperature. The difference 
in the relative offset frequency between the upshifted and downshifted lines allows a 
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determination of the bulk electron drift velocity. This velocity can be compared with 
the bulk velocity estimated from the ion line measurement to calculate the relative 
drift between the electrons and ions (i.e., current). 

In summary, we have presented the characteristics of the normal incoherent scat- 
ter power spectrum. The spectral features of the recently observed enhanced radar 
backscatter, examined in the next section, exhibit significant departures from the 
normal spectrum. 


5.2 Characteristics of Enhanced Radar Backscat- 
ter (ERB) 


Intense coherent radar returns from the topside of the ionospheric F region have 
been a regular feature of the Millstone Hill UHF incoherent scatter radar data for 
many years. Because these returns are characterized by the large radar cross sections 
and narrow spectral widths which can result from satellite penetration of the 1° radar 
beam, they have, until recently, been dismissed solely as hard target contamination of 
the incoherent scatter data. Foster et. al [1987] noted an anomalously high occurrence 
of intense radar echoes in the vicinity of the mid-latitude ionospheric trough. The 
unique geophysical conditions associated with the trough, such as field-aligned current 
activity and steep plasma density gradients, suggest that some of the coherent returns 
may be caused by enhanced ionospheric plasma density fluctuations, rather than 
orbiting objects intersecting the radar beam [FOS88]. 

The generation of such spectra may be attributed to current- driven ionospheric 
plasma processes [FOS88]. The processes described by Kindel and Kennel [1971] 
and Rosenbluth and Rostoker [1962], for example, predict the amplification of ion 
acoustic waves in the presence of intense currents and explain several features of the 
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observed enhancements, as discussed in Chapter 7. Here we present evidence derived 
from several sources for the observation of stimulated ion acoustic wave growth with 
the Millstone Hill UHF radar. These include statistical analyses of large data sets, 
multi-diagnostic measurements, and the observation of ERB in zenith experiments. 

5.2.1 Satellite Contamination and Statistical Analyses of 
ERB 

The task of separating satellite returns from true ionospheric coherent echoes is diffi- 
cult because the integrated spectral features of both axe essentially identical. This is 
illustrated in Figure 5-2, where the power spectra of the normal incoherent scatter ion 
line, an intense return from an enhanced ion line, and spherical satellite number 9636 
are plotted on the same frequency scale for comparison; each plot is self-normalized. 
Even on a pulse-to-pulse basis satellite returns may vary by as much as 20 dB due to 
reflecting surface irregularities and rotation effects [5. Sridharan , private communi- 
cation, 1988]; the size of this variation is comparable to the magnitude of the largest 
geophysical ERB events. Thus, the similarity between geophysical ERB spectra and 
satellite spectra makes unambiguous classification of individual ERB events difficult. 
Severed approaches to the resolution of this problem have been employed; the most 
successful involves the statistical reduction of large data sets, as described below. 


The Satellite Catalog 

A large number of satellites pass through the Millstone I.S. Radar field of view during 
any particular experiment; their spectral signatures constitute the primary source of 
coherent echo contamination in the topside ionosphere. At Millstone Hill we are 
fortunate to have access to a complete catalog of the known orbiting space objects. 
An altitude distribution of orbiting objects generated from this catalog is shown in 
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Figure 5-2: UHF Power spectrum of a) Normal incoherent scatter ion line; b) En- 
hanced ion line; and c) Spherical satellite ^9636. The amplitude on each plot is 
self-normalized. 
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Figure 5-3: Altitude distribution of satellites in the Millstone Hill radar field-of-view 
for a 6° elevation angle. 

Figure 5-3. 

The majority of the satellite population exists at altitudes above about 500 km; 
layers of debris can be identified near 600 km and 800 km, with a well defined max- 
imum near 1000 km. The distribution underscores the problem of removing hard 
target contamination: The Millstone Radar cannot diagnose the active auroral re- 
gions north of 60° latitude at altitudes less than 500 km. Thus we are forced to 
sample the active region of interest at altitudes where the greatest number of con- 
taminants are located. A number of techniques utilizing the satellite catalog have 
been employed to distinguish satellite echoes from those of geophysical origin. 

Specific Events One way to employ the catalog is to simply check specific records 
of the radar data containing coherent peaks suspected to be geophysical ERB returns. 
The time and location of a peak’s occurrence can be checked against the catalog to 
determine if a satellite with the appropriate Doppler velocity was present. In applying 
this method one must be careful to properly define the uncertainty in the space-time 
window which must be specified for the coherent echo s location. 
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A single Millstone data record is typically constructed of power measurements 
which are integrated for 15-30 seconds, during which time the radar may scan up to 
five degrees in azimuth. Spatial range smearing due to the radar pulse length adds 
up to 150 km uncertainty in range. Additional error is introduced by the finite radar 
beamwidth; the large radar cross section of many satellites requires one to define 
the “effective beamwidth” as including the sidelobe pattern well beyond the common 
definition of the half-power width. For the Millstone Hill steerable UHF antenna, this 
means expanding the traditionally defined 1° beamwidth to some 2°-3°. The largest 
uncertainties in the azimuth, range, and time of the space-time window specified for 
a given echo are then, respectively, ± 5°, 150 km, and 30 seconds, enough time for a 
satellite to travel some 250 kilometers. 

In reality, of course, the radar beam occupies only a small fraction of the total 
window at a given moment in time; to insure a fair check, however, the entire window 
must be specified for the satellite catalog to check a particular data record. The result 
is that the catalog often predicts the penetration of the relatively large space window 
by several hard targets within the 30 second record interval, even when no evidence 
of a corresponding coherent peak is seen in the data. 

If the catalog predicts that no suitable satellites were present, one can at least be 
confident that the assumption of the peak’s ionospheric origin has not been ruled out. 
While this use of the satellite catalog is perhaps the most obvious, the technique can- 
not provide definitive proof that a given radar echo was generated by the ionosphere 
because the satellite catalog is not 100% reliable. 

The catalog is not completely reliable for several reasons. The large number of 
objects cataloged (over 7000 presently) makes it very difficult to provide frequent 
element set updates on all objects. The low altitude orbits occupied by the satel- 
lites penetrating the Millstone Hill radar beam degrade fairly rapidly, so that catalog 


83 


projections of the locations of some objects may be significantly inaccurate. Fur- 
thermore, there are numerous classified satellites not included in the catalog used for 
these studies. Finally, there is a population of uncataloged small orbiting debris. 

Statistical Studies The uncertainties involved in applying the satellite catalog to 
check specific events limit the usefulness of that technique, as noted above. The 
factors that render the catalog fallible for selected records, however, become unim- 
portant if one applies the catalog in a temporally statistical manner. A calculation 
of the distribution of low altitude targets < 110 minutes) passing through a 

defined space window over a period of 24 hours or so, will not depend significantly 
on the accuracy of a particular satellite’s ephemeris data. Assuming they represent 
a small percentage of the total satellite population and exhibit a similar altitude dis- 
tribution, classified objects no longer constitute a serious source of error, either. The 
contamination due to the uncataloged debris must be examined in more detail. 

Only about five percent of the over 7000 space objects cataloged are active or 
working satellites. About 50% of the known objects orbiting the earth consist of 
fragments resulting from the breakup of artificial satellites. The primary causes of 
satellite breakups are propulsion related malfunctions and intentional detonation, 
although collisions are suspected for unexplained breakups as well [ Johnson , 1985]. 
It is believed that the uncataloged population of fragments from these breakups may 
be up to 2.5 times greater than the total cataloged population [ Kessler , 1985]. 

Statistical studies, however, will not be modified significantly by the presence of 
the orbiting debris for two reasons. The first is that the uncataloged population con- 
sists of fragments smaller than about 10 centimeters diameter. Objects of this size 
in the topside of the ionosphere are too small to cause greatly enhanced coherent 
echoes, particularly for radar wavelengths much greater than the fragment diameter. 
The second reason is that the distribution of fragments generated by a breakup re- 
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mains centered in height at the original orbit altitude. Thus the shape of the altitude 
distribution changes very little; only the number of objects varies. 

Doppler Shift Discrimination The catalog has been employed to calculate the 
passage of satellites through an azimuth window which is left “open” for a sufficient 
period of time to allow several orbital revolutions, typically 24 hours. In Figure 5-4a a 
histogram of the number of cataloged objects is plotted as a function of their line-of- 
sight velocity to the Millstone radar for a low elevation azimuth scan between 342°-27° 
over a 24 hour period. The spikes located at about ±6 km/sec are due to the large 
population of satellites in so-called “polar” orbits, high inclination orbits between 
800-1000 kilometers altitude. The radar measures a large line-of-sight velocity as it 
scans to the north and detects these objects as they pass over the polar cap. The 
velocity measured by the radar is actually determined from the Doppler shift of the 
scattered signal received by the antenna. For the backscatter geometry employed 
here, the shift in kHz is given by Af = 2v/X, where v is the satellite velocity in 
km/sec and A the radar wavelength in meters. A simple conversion for the 68 cm 
wavelength used by Millstone Hill is Af — 3u. 

Figure 5-4b shows the distribution of coherent peaks observed at Millstone Hill 
for an experiment run 13-15 January, 1988 scanning through the same azimuths as 
the simulation 5-4a. A comparison of the two nearly identical distributions suggests 
that many, if not all, of the peaks observed during this experiment were caused by 
satellites penetrating the scanning radar beam. Magnetic conditions during 48 of 
the 72 hours of data analyzed were very disturbed; Kp values of 7 or higher were 
registered for a 9 hour period, and remained at or above 6 for 15 consecutive hours 
during the middle period of the experiment. 

The consideration of magnetic activity is important because of its correlation with 
field-aligned currents. The Kp indices referred to here are global classifiers of the am- 
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plitudes of magnetic variations derived every three hours based on observations from 
a number of stations located over a wide range of latitudes. The quasi-logarithmic 
scale varies from 0 to 9 in 27 increments; the scale is skewed towards the low end of 
these values in that magnetic fluctuations corresponding to 3 or less occur 90% of the 
time, whereas fluctuations of Kp=6 or greater represent less than 1% of the total dis- 
tribution [MayaucL, 1980]. Given the sustained periods of very high activity occurring 
during the experiment of JAN88, we might expect a correspondingly high intensity 
of field-aligned currents, producing significant numbers of ERB events. Reasons for 
the apparent lack of observed ERB in this experiment will be discussed shortly. 

A second experimented simulation and corresponding observations are presented 
in Figure 5-5. As in Figure 5-4, the signature of the polar orbiters at ±6 km/ sec is ev- 
ident in both 5-5a and b. However, in Figure 5-5 an even larger number of peaks were 
observed between ±1-2 km/sec, the ion sound speed regime in topside ionospheric 
plasmas. This remarkable feature has not been observed in velocity distributions 
constructed from other Millstone Hill radar experiments. Clearly the satellite catalog 
does not predict the large population of peaks observed at these Doppler velocities. 
If this data set represents a statistical average, we can feel confident that a portion 
of the observed ion acoustic peaks were generated by ionospheric processes. 

The satellite simulations shown in Figures 5-4a and 5-5a were generated by contin- 
uously monitoring all azimuths within a defined window simultaneously, in so-called 
open window simulations. Because the simulated coverage area is large, 12 hours 
(t > 67-OT.bit) is a sufficient time to run the open window simulation. Simulating longer 
periods of time does not affect the shape of the calculated distribution of detected 
hard targets; it only increases the number of “hits” in the radar beam. 

In an actual experiment, of course, the beam scans through the window and 
monitors only a small fraction of the total azimuth extent at any moment in time. The 
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Figure 5-4: Distribution of Doppler velocities for a) satellites penetrating the radar 
beam between 342°-27° AZ, 4° EL, and b) Actual observations from 13-15 Jan. 1988 
for the same pointing angles. 
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Figure 5-5: Distribution of Doppler velocities for a) satellites penetrating the radar 
beam between 270° -90° AZ, 6° EL, and b) Actual observations from 8-10 November 
1987 for the same pointing angles. The large number of peaks observed at ion acoustic 
velocities is not predicted. 
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Figure 5-6: Graphical illustration of procedures used for a) satellite simulations where 
the entire space window is continuously monitored, and b) actual data acquisition via 
scanning radar beam in experiment. 

difference between the simulations shown thus far and actual experiments is illustrated 
graphically in Figure 5-6. Hence, the number of satellites actually detected in a 12 
hour observational period is much smaller than the number “seen” by the catalog for 
the same 12 hours, and a comparison of the corresponding observed distribution may 
not be statistically valid. 

The minimum time required for a valid statistical comparison of observed and 
simulated coherent peak velocity distributions was determined by simulating a com- 
puter radar beam scanning through a dynamic satellite field. In the simulations, the 
table of satellites passing through the azimuth window were sorted according to time, 
and a radar scan was initiated at a given time and starting position to sample the 
satellite data base. The simulated scan consisted of a 2° x 2° wide beam which was 
swept 10°/min in azimuth, the same azimuth scan rate used in most experiments at 
Millstone Hill. This process was repeated several times using different initialization 
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parameters for periods of 12, 24, and 48 hours. 

It took 4-5 times more scanning hours to acquire as many targets as the corre- 
sponding open window simulation. The 12 hour simulations produced distributions 
which were usually similar to the open window distribution, but in some cases rela- 
tively large disparities were evident. The 24 hour simulations improved results. None 
of the 24 hour distributions, for example, predicted more peaks at ion acoustic ve- 
locities than at polar orbiter velocities, as was observed in the NOV87 experiment. 
Nevertheless, deviations from the open window distribution were large enough to 
discourage quantitative comparison. After 48 hours, however, the scan simulations 
matched the shape of the open window results nicely, regardless of starting time or 
position within the satellite data base; increasing the time further simply produced 
more events without noticeable distortion of the distribution. 

The data sets used to produce the observed distributions presented in Figures 
5-4b and 5-5b contain at least 60 continuous hours of radar observations, and should 
therefore satisfy the statistics of the satellite catalog simulations. Comparing the 
simulations with the observed distributions we conclude that a number of the co- 
herent peaks observed at ion acoustic velocities between 1—2 km/sec in the NOV87 
experiment result from geophysically induced ERB. The levels of global magnetic 
activity during the period varied from quiet to unsettled to moderately disturbed 
(0 <Kp< 5"). Extrapolating the number of ERB events occurring during this mod- 
erate activity to the extremely disturbed periods of the JAN88 experiment, we would 
expect to find much stronger evidence for geophysical ERB in the January velocity 
distribution. The lack of such evidence can be explained by considering the results 
of a statistical analysis of a third experiment. 

Altitude Discrimination Coherent echo distributions were constructed from ap- 
proximately 80 hours of data acquired during an experiment on 6-10 March 89, which 
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was known to contain magnetically disturbed periods. The velocity distribution of the 
enhanced echoes observed during the experiment, however, very nearly matched the 
corresponding distribution generated by the satellite catalog. The observed altitude 
distribution of the enhanced returns, shown in Figure 5-7, also looks quite similar 
to the simulated satellite distribution shown in Figure 5-3, except at the lowest al- 
titudes, where a surprisingly large number of peaks are observed in range gate 10, 
corresponding to 330 km altitude. 

The number of orbiting objects at altitudes of about 300 km or less represents less 
than .5% of the total satellite population. In the observed distribution, however, such 
events account for nearly 1.5%. For this experiment that translates into precisely 18 
events, a small number which is perhaps not statistically robust, despite 80+ hours 
of observations. The significance of the low altitude anomaly increases when one 
recognizes that the number of events near 300 km is comparable to, or greater than 
the number of peaks detected at altitudes of 500 and 540 kilometers, respectively. An 
examination of the spectral characteristics of each low altitude enhancement provided 
further evidence of a geophysical source for some of the enhancements. 

Satellites at low altitudes must travel faster than orbiting objects at higher al- 
titudes; a satellite simulation for the 6—10 March azimuth scans shows that about 
85% of the time, the line-of-sight Doppler velocity detected by the radar from these 
objects exceeds 4 km/sec, and is often higher than 7 km/sec. In the March experi- 
ment, however, 5 of the 18 low altitude enhanced peaks exhibited Doppler velocities 
at the local ion acoustic speed, usually between 1 and 2 km/sec in the F region. This 
is almost double the expected percentage of ion acoustic enhancements, and seems 
unlikely to be a random occurrence caused by hard targets. The occurrence of these 
events with regard to magnetic activity and other geophysical considerations is also 
indicative of ion acoustic enhancements through natural processes. 
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ALTITUDE DISTRIBUTION OF OBSERVED ERB: 6-10MAR89 



Figure 5-7: The altitude distribution of observed coherent echoes from the Millstone 
6-10 March 1989 experiment. Altitude is scaled by range gates, which are separated 
by 40-50 kilometers. Some gate numbers followed by their corresponding altitude 
are: 10/330, 15/500, 20/700, and 25/975. 

This March experiment time period was very interesting magnetically, being more 
disturbed than the 8-10 NOV87 period, but significantly less disturbed than the mon- 
ster storm levels recorded during 13-15 JAN88. The hourly AE index is plotted for 
the duration of the experiment in Figure 5-8. The AE index provides a useful mea- 
sure of auroral magnetic activity. The figure shows that March 6—7 were moderately 
disturbed before activity decreased steadily through the latter part of the 7th and 
March 8. At 18:00 UT on the 8th a major storm commenced; disturbance levels de- 
creased soon after the major impulse, but remained active throughout most of March 
9 when another surge of activity occurred, lasting just til 24:00 UT. 

The height of the maximum density in the F region, hmF2, is closely related to 
the magnetic activity (see Figure 5-9). Under quiet conditions, hmF2 exhibits a reg- 
ular diurnal variation, descending during daylight hours as photoionization builds up 
the plasma density at lower altitudes, and ascending again after dusk when recombi- 
nation processes neutralize the low altitude plasma more rapidly. A typical variation 
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Figure 5-8: Values of magnetic index AE during the magnetic disturbances of 6-10 
March, 1989. Note the sudden storm commencement at 18:00 UT on March 8 [after 
Buonsanto et al ., 1990]. 

of this nature is evident during the period of decreasing magnetic activity spanning 
March 7-8. Under magnetically disturbed conditions hmF2 serves as an excellent 
indicator of the passage of the mid-latitude electron density trough and the equa- 
torward expansion of the auroral oval [Buonsanto , et al., 1990]. In the trough itself, 
the hmF2 value increases sharply, as occurs in the early morning hours of March 7, 
March 9, and March 10. 

The low altitude ion acoustic enhancements occurred at 7MAR-1:44 UT, 9MAR- 
3:34 UT, 9MAR- 11:52 UT and 10MAR-3:00 UT. Their positions are indicated on 
Figure 5-10, which is a plot showing the variation of the peak ionospheric density, 
NMF2, as a function of time and latitude. All the enhanced peaks occur either in 
the ionospheric trough or on the trough’s equatorward edge in a region of strong 
density gradients associated with field-aligned Birkeland currents [Ungstrup, et al., 
1986]. Except for the event at 11:52 UT on March 9, all the events occur just prior 
to the arrival of the trough over Millstone Hill. The event at 11:52 UT is actually 
occurring under similar geophysical conditions, except that it is located on the edge 
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Figure 5-9: The height of maximum electron density in the ionosphere for the 6-10 
March period. The sharp increases in altitude correspond to the approach of the 
mid-latitude trough [after Buonsanto et al . , 1990]. 

of a receding, rather than approaching, density depletion. 

The approach of the trough signals the expansion of the auroral oval and the 
penetration of large electric fields to lower latitudes. These fields are responsible 
for the very large westward E x B plasma drifts observed in the early morning on 
March 9. The enhanced low altitude ion line detected at 3:34 UT on March 9 is 
shown in Figure 5-11; the large shifts evident in the normal IS spectra at higher 
altitudes correspond to westward drifts of up to 1.3 km/sec. All-sky images revealed 
that prominent red aurora were present north of the Millstone Hill radar at this 
time, in the same location where the ERB occurred; later in the morning a stable 
red auroral (SAR) arc persisted to the north for about two hours [M. Buonsanto , 
private communication, 1991]. These geophysical signatures verify the presence of 
auroral features at latitudes as low as 44° on March 9. Similar features, though not 
as intense, were present during the early morning hours of March 7 and 10 as well. 
By contrast, the post-dusk period on March 7/8 did not display a well-defined trough 
or auroral characteristics. The density depletion evident after midnight on this day 
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Figure 5-10: Variation of peak electron density with latitude. Low altitude ERB 
events are indicated with white dots. Low density regions are shaded black [after 
Buonsanto et al., 1990]. 
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Figure 5-11: Ion line enhancement observed with the Millstone Hill UHF radar at 
03:34 UT on 9 March 89. At higher altitudes the shifted spectra give evidence of 
large westward E x B drifts. 

is associated only with the normal diurnal decrease occurring after sunset. No low 
altitude ion acoustic ERB events were recorded during this time period. 

We have established the presence of active auroral features during the periods 
when low altitude ERB events were observed in the 6-10 March experiment. These 
events are attributed to the enhancements of ion acoustic waves in the presence of 
naturally occurring intense field-aligned currents. Normally the Millstone Hill Radar 
can diagnose these remote high latitude active regions only at altitudes above 500 km, 
where the large satellite population makes unambiguous identification of ERB events 
difficult. During the magnetic disturbances of March 6/7, March 8/9 and March 10, 
the equatorward expansion of the auroral oval carried it to within a few degrees of 
Millstone’s latitude. This made it possible for low elevation- angle azimuth scans to 
sample the active region at much closer ranges, and correspondingly lower altitudes, 
than normal. The low altitude ERB events were then relatively easy to separate 
from hard targets because so few satellites occupy these very low altitude orbits. The 
variability of the ionosphere played an important role in interpreting the results of 
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the distribution studies for this experiment. 

The level of variability in the 13-15 JAN88 experiment was about an order of 
magnitude greater than that experienced during 6-10 MAR89. By extending the 
auroral dynamics arguments used to explain the results of the MAR89 experiment, 
we can arrive at a rational explanation for the total lack of statistical evidence of 
geophysical ERB during the JAN88 experiment. In MAR89, ionospheric variations 
were characterized by brief periods of intense magnetic activity resulting in the equa- 
torward expansion of the the auroral oval to latitudes slightly north of Millstone Hill. 
The periods of intense activity were then followed by somewhat longer periods of 
decreased activity, providing a relaxation period for the ionosphere to nearly return 
to its normal structure. 

After the first day of the 13-15 JAN88 experiment, however, periods of extreme 
activity were sustained for many hours. The auroral region, containing the intense 
current activity critical for ERB excitation, actually passed to the south of Millstone 
Hill and did not recover to its usual latitudes during the duration of the experiment. 
The radar scan cycle employed, a 4° elevation-angle, limited azimuth scan to the 
north, was not advantageous to the detection of ERB under these circumstances. 
Zenith measurements or even scans to the south of Millstone are appropriate for such 
studies under very disturbed conditions. In conclusion, ERB events were not detected 
in the JAN88 experiment because the auroral zone passed to the south of Millstone 
Hill due to extreme levels of magnetic activity, while the radar continued to scan at 
low elevation angles to the fax north. 

Finally, a large data set acquired on 6-10 October 86 was analyzed by the same 
statistical technique. The observed velocity distribution closely resembled the satellite 
distribution, similar to the comparison shown in Figure 5-4. Furthermore, an analysis 
of the altitude distributions yielded no evidence of anomalous ERB events. The global 
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Kp index indicated that the period was one of the most quiet imaginable, with 60 of 
the 87 observing hours having Kp of 1 + or less and a maximum Kp of 3' occurring for 
one 3-hour period. Under these conditions, intense current activity is not expected at 
any latitude within Millstone Hill’s field-of-view. Consistent with that interpretation, 
we conclude that the lack of evidence of ERB events in the statistical analyses for 
this experiment indicates simply that no appreciable ERB events occurred during this 
time period because of the very low levels of magnetic activity. 

5.2.2 MICAD Experiments 

Beginning in April 1989, a series of multi-diagnostic experiments, known as MICAD . 
involving the Millstone Hill UHF radar, a subset of the CA NOPUS system of in- 
struments (BARS radar, magnetometer array, meridian photometers) and the DE-l 
satellite was conducted to provide additional information about ion acoustic radar 
backscatter enhancements due to ionospheric processes. In these experiments, as 
illustrated in Figure 5-12, the Millstone Hill radar was used to track the geomag- 
netic field line above the DE-l satellite as it flew through perigee within the radar’s 
field of view; CANOPUS instruments were utilized as diagnostics of the background 
ionosphere. Ultimately, the goal was to simultaneously detect the signature of large, 
field- aligned ionospheric currents with remote sensors (radars, etc.) and in situ in- 
struments aboard DE-l. 

Beginning in April 1989, a total of seven MICAD overflight experiments were 
conducted on 4/9, 5/2, 5/22, 6/5, 6/19, 8/8, and 8/29. During the MICAD exper- 
iment on 22 May 1989, the Millstone Radar detected a 7 dB enhanced backscatter 
return at the ion acoustic frequency at 23:57:46 UT; the power enhancement at 1000 
km is shown in Figure 5-13. The power plot shows a distinct F region trough at 
the latitude of the ERB event. In the scan shown in the figure the radar moved 
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Figure 5-12: Conceptual design of MIC AD experiments. Remote diagnostics and in 
situ instruments were to obtain simultaneous measurements on the same magnetic 
field line. 
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from AZ=289° at 23:54:52 UT to AZ=343° at 23:59:08 UT. The enhancement oc- 
curred at an azimuth of 332°. The DE-1 satellite passed through the same region 30 
seconds later and recorded an in situ magnetometer fluctuation corresponding to a 
field-aligned current density > 50 /xamps/m 2 [J. Slavin , private communicaton, 1989]. 
The magnetometer trace and the enhanced UHF radar spectra are shown in Figure 
5-14. The decrease in the magnitude of DBPHI evident just after 23:57:00 UT in the 
figure indicates the satellite’s entry into the Region II field-aligned current system. 
The boundary between the Region I currents and the Region II currents is reached 
at about 23:58:10 UT, when the ~ 50 /xamps/m 2 current spike was detected. The 
onboard magnetometer continued to exhibit large fluctuations until about 23:59:00 
UT; two of the additional fluctuations at 23:58:20 UT and 23:58:40 UT correspond to 
current densities in excess of 25 /iamps/m 2 . Because the radar was actually tracking 
about 30 seconds ahead of the satellite’s position, the measurements cannot be clas- 
sified as truly simultaneous. However, this experiment gives firm evidence of intense 
field-aligned currents associated with the occurrence of ERB near the ionospheric 
trough as was suggested by FOS88. 

An examination of data from two consecutive passes of the DMSP-F8 satellite 
(the first 42 minutes prior to, and the second 56 minutes after the DE-1 overflight 
in the same longitude band) reveal that the region producing the enhanced radar 
backscatter was characterized by localized intense, soft (E < 1 keV) electron pre- 
cipitation features. Comparison of the measured electron and ion fluxes show that 
current densities in excess of 60 /iamps/m 2 were detected on the first pass [W. Denig, 
private communication, 1989]. Data from this pass formed the basis for deriving the 
suprathermal contribution to the electron distribution function utilized in numerical 
calculations in Chapter 7 (see Figure 7-7). 

Particle flux spectrograms from a northern polar pass of the NOAA-10 satellite 
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Figure 5-13: Power enhancement observed during May 22 MICAD experiment. The 7 
dB enhancement is located at 1000 km altitude on the edge of the ionospheric trough. 
A two msec radar pulse was used in the scan. 


exhibit a sharp local maximum in the low energy electron flux at the same latitude 
and only 8 minutes prior to the DE-1 satellite and Millstone Hill radar observations. 
A map illustrating the geographic and temporal relationship of the various measure- 
ments is shown in Figure 5-15. These in situ data establish the presence of discrete 
low energy particle precipitation events carrying large current densities in the region 
of enhanced UHF radar backscatter. 

It should be noted that limited particle flux information is available from the 
DE-1 satellite. The electron flux instrument onboard is no longer functional. The 
energy spectrum of ions was monitored, but no signature was evident in association 
with the large magnetometer variations [ W. K. Peterson , private communication, 
1989]. This is consistent with the DMSP data, where the energy flux of the electron 
population increases dramatically while the ion flux remains essentially unchanged. 
Both measurements indicate that the currents are being carried by suprathermal 
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Figure 5-14: Power spectra of ion acoustic ERB observed with the Millstone Hill UHF 
radar and the corresponding magnetometer trace of the DE-1 satellite. Magnetic 
fluctuations (DBPHI) are in units of nanotesla. 
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Figure 5-15: Map showing geographic relationship of multiple satellite measurements 
of intense field-aligned current activity, UHF ion acoustic ERB, and the Millstone 
Hill radar. 


electrons rather than ions. 

Measurements from the CANOPUS diagnostic group of instruments were of lim- 
ited application in this particular MICAD experiment. The observations were outside 
the field-of-view of the BARS radar system, and the local time was too early to ac- 
quire all sky camera images of visible airglow or photometer data. The meridional 
magnetometer chain (MARIA) approximately 320 km west of the UHF ERB observa- 
tion was operating, and the magnetometer traces from the nearest stations are shown 
in Figure 5-16. The ground stations show a reasonable level of activity (Kp=4) prior 
to the Millstone/ DE-1 observations, but conditions are relatively quiet during the 
actual overflight (about 2 hours after start time). Global magnetic activity usually 
serves as an indicator of intense field-aligned currents, but it does not constitute a 
prerequisite for their occurrence. Johnstone and Winningham [1982] reported several 
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satellite observations of suprathermal electron bursts carrying field-aligned currents 
in the auroral zone during very quiet magnetic conditions. 

The satellite catalog was consulted for the presence of hard targets in the scattering 
volume during the 30-second integration period in which the ion acoustic ERB event 
was recorded. A generous scattering volume was assumed, using the actual elevation 
and azimuth angles of the radar ±2°, to minimize the uncertainties due to errors 
in the catalog predicts. Three candidates were found. Of these three, two were 
expected to produce large line-of-sight Doppler shifts of about 6 km/sec. The other 
possible candidate's Doppler velocity was in the opposite direction of the observed 
phase velocity of the enhanced ion acoustic peak. Clearly, none of the predicted 
hard targets could have caused the asymmetric UHF power spectrum observed in 
this MIC AD experiment. 

5.2.3 Observation of ERB in Zenith Experiments 

As has been established, the contamination due to satellites in the investigation of 
ERB is a serious problem in general. However, hard target signatures can be filtered 
out effectively by pointing the radar directly overhead. In the zenith position the 
beam is nearly transverse to the direction of motion of satellites penetrating the 
scattering volume, and their associated Doppler shifts are thus restricted to nearly 
zero magnitude. Satellite catalog simulations based on a 2° x 2 wide beam show 
that the maximum observable line-of-sight velocity possible in such experiments does 
not exceed 250 m/sec. Zenith observations of spectral peaks with greater Doppler 
shifts (e.g., F - region ion acoustic speed: C a ~l-2 km/ sec) can then be attributed to 
geophysical sources. 

Figure 5-17 shows an observation of ion acoustic wave growth acquired with the 
Millstone Hill 67-meter fixed position zenith antenna. The Doppler shift of the peak 
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Figure 5-16: Magnetic variations recorded by the CANOPUS chain of magnetometers 
during the May 22, 1989 MICAD experiment. The DE-1 and Millstone Hill obser- 
vations of ERB were recorded at 23:58 UT, during the relatively quiet period nearly 
two hours after the start time on the plots. 
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corresponds to a 1 km/sec velocity, about four times greater than the maximum 
observable shift from actual satellites. In addition, the enhanced return is distributed 
in altitude to a greater extent than is possible for an orbiting object given the 150 km 
pulse smearing along the radar beam. This zenith observation, reported by FOS88, 
was taken as proof that a portion of the Millstone ERB returns cannot be caused by 
orbiting objects. 

The scan cycle for the radar in this experiment consisted of 4 consecutive zenith 
measurements followed by a complete elevation scan, 2 more zenith records, and a 
180° low elevation-angle azimuth scan. Only six minutes out of every hour are spent 
in the zenith position, and only half that time is used to diagnose the high altitude 
region above 700 km. Unfortunately, the spectra shown in Figure 5-17 represent 
the last zenith measurement preceding the beginning of an elevation scan lasting six 
minutes. When the radar returned to the zenith position, the echoes were gone. 

Observations like this at Millstone Hill are rare for two reasons. As was noted 
previously, the enhancements are thought to be associated with intense field-aligned 
currents. While such events are not extraordinary at auroral latitudes, they occur 
at mid-latitudes very infrequently, approximately 10—15 days a year based on the Kp 
magnetic index during the event shown here [ MayaucL, 1980]. On those days when the 
general levels of disturbance are sufficient, the phenomenon directly over the radar 
may last for only minutes or seconds and be missed entirely [Rietveld et ai, 1991], as 
was nearly the case in this experiment. 

Typically only a small amount of Millstone’s total data acquisition time is used 
for zenith observations. Few dedicated zenith experiments have been conducted at 
Millstone Hill in the past; two separate 24- hour experiments run recently yielded no 
significant results. The zenith data obtained in most previous experiments consists 
of a few local measurements made once or twice per hour in between azimuth and/or 
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Figure 5-17: Millstone Hill zenith observation of ion acoustic ERB. The frequency 
offset corresponds to a velocity of 1 km/sec. [after FOS88]. 
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elevation scans. Considering the low occurrence frequency of the geophysical event 
and the sparse number of zenith observations, the probability of detecting ion acoustic 
enhancements directly above Millstone Hill is small; zenith experiments conducted 
solely for that purpose may not provide an effective approach for investigating ERB 
at mid-latitudes. 


5.3 Summary 

An attempt has been made to provide conclusive evidence of geophysically-produced 
enhanced ion acoustic radar backscatter observed at Millstone Hill. 

Because of the spectral similarities of satellites and enhanced ion acoustic waves, 
it is difficult to distinguish individual measurements unambiguously. Satellite sim- 
ulations have been employed to make statistical comparisons with large data sets. 
Results from four types of experiments have been presented. Under extremely dis- 
turbed and extremely quiet conditions no ERB events were detected by the statistical 
analysis. In the quiet case there are most likely no occurrences of ERB. In the very 
disturbed case, however, pointing the radar to the north proved ineffective because of 
the auroral region’s passage to the south of Millstone. ERB events were found at low 
altitudes during a lesser storm period when the auroral oval expanded to within a few 
degrees north of Millstone. These events have been attributed to the current- driven 
enhancements of ion acoustic waves in the temporarily nearby auroral zone. Finally, 
a moderately disturbed period was examined. The velocity distribution exhibited 
many more peaks at ion acoustic frequencies than predicted by the satellite catalog. 
These peaks are believed to represent geophysical ERB events occurring at relatively 
high latitudes in the usual auroral zone. 

Through a series of multiple diagnostic experiments, we have established the asso- 
ciation of enhanced backscatter events with intense field-aligned currents. In one case, 
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a nearly simultaneous measurement of UHF ERB and a field-aligned current carried 
by suprathermal electrons were obtained. Data from other satellites confirmed that 
the region contained large fluxes of low energy electrons, and that such conditions 
persisted for tens of minutes. These soft precipitating electrons are believed to play 
an important role in the generation of ion acoustic ERB. 

Finally, the utility of zenith-pointing experiments was exploited to filter out hard 
target contamination while acquiring ion acoustic ERB. A single event has been iden- 
tified with a Doppler shift of 1 km/sec, the local ion sound speed in the plasma. 
Satellites passing through the beam exhibit Doppler shifts of less than 250 m/sec. 
These measurements are rare at mid-latitudes because intense field-aligned currents 
occur there very infrequently (perhaps briefly on 10—15 days/year during solar max- 
imum), and the radar is not usually pointed in the zenith direction while acquiring 
data. Results from two 24-hour zenith experiments yielded no evidence of ion acoustic 
ERB. 
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Chapter 6 


ERB Observations from Other 
Radar Sites 


The existence of geophysically induced enhanced coherent radar backscatter from the 
topside ionosphere was first reported by FOS88 based on a study of “hard target 
like returns observed with the Millstone Hill UHF Radar as described in Chapter 5. 
The enhancements are believed to be associated with intense field-aligned currents 
in the high latitude ionosphere. This high-latitude region is normally accessible to 
the Millstone Hill radar beam at very large aspect angles relative to the geomag- 
netic field, effectively limiting the range extent over which the field-aligned echoes 
may be observed and leading to ambiguities (with actual hard targets) in the data 
interpretation for individual coherent echo events. The locations of the EISCAT and 
Sondrestrom radars (see Table 4.1) suggest that more favorable ERB observing con- 
ditions may exist at those sites. Presented here is a brief overview of significant 
observations and experiments involving these radars which confirm and extend the 
known characteristics of ERB first recognized at Millstone Hill. 
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6.1 Observations at EISCAT 


The EISCAT (European Incoherent SCATter) Scientific Association operates a tri- 
static UHF radar system consisting of a 933 MHz 2.5 MW radar located at Troms0, 
Norway (69.7°N lat, 19.2°E.long) and two remote receivers at Kiruna, Sweden (67.8°, 
20.4°) and Sodankyla, Finland (67.4°, 26.6°) . The three sites utilize identical parabolic 
antennas 32 meters in diameter to simultaneously determine full three-dimensional 
vector velocities at selected heights in the ionosphere. In addition to the UHF radar, 
a 224 MHz radar using a 40 x 120 m parabolic cylinder antenna is located at Troms0, 
providing multi-wavelength measurement capability locally. 

6.1.1 Initial Statistical Approach 

The investigation of ERB at EISCAT was first undertaken by Schlegel and Moorcroft 
[1989] who applied a statistical approach, similar to that used at Millstone Hill, 
to 73 hours of data acquired with the radar pointed up the local magnetic field 
line (AZ=182°, EL=77.5°). EISCAT regularly runs so called Common Programs for 
data acquisition which employ this radar look angle for geophysical reasons, such 
data constitute a majority of the available data base. Attempting to discriminate 
geophysical events from hard targets, Schlegel and Moorcroft’s analysis identified 
and correlated 209 enhanced events with various combinations of occurrence time, 
altitude, velocity, intensity, and spectral characteristics. Achieving ambiguous results, 
they concluded that the ERB observed at EISCAT was probably due to a combination 
of satellites and system effects not clearly understood. 

The satellite explanation was supported in part by a comparison of the distribu- 
tion of observed peak velocities with the distribution predicted for hard targets by a 
satellite catalog. While this comparison was employed successfully with the Millstone 
Hill data, it fails to discriminate ERB events in this case because the field-aligned 
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pointing direction at EISCAT produces radar line-of-sight Doppler velocities from 
the large population of polar orbiting satellites in the range of typical ion acoustic 
velocities between 1-2 km/sec (see Figure 6-1). 

However, a marked difference in the respective height distributions of satellites 
and the observed peaks was discovered. Shown in Figure 6-2 is the unexpectedly 
large population of observed low altitude events. Aware that reports of Millstone Hill 
ERB were restricted to the topside ionosphere above 400 km [FOS88] and lacking 
additional evidence of geophysical processes, Schlegel and Moorcroft were unable to 
resolve the disparity in the altitude distributions, leaving the question open for further 
interpretation. Recently, data reported by Rietveld et al. [1991] (henceforth RD91) 
unambiguously verified the detection of geophysically induced ERB with the EISCAT 
radar system. 

6.1.2 Range Extended UHF Echoes 

Like their counterparts observed at Millstone Hill, the geophysical ERB echoes de- 
tected at EISCAT were also characterized by large cross section enhancements (up 
to 20 dB) and narrow spectral widths similar to those caused by hard targets. These 
grossly asymmetric spectra cannot be fitted to derive ionospheric parameters by the 
standard incoherent scatter analysis programs. When the radar is pointed up the 
magnetic field line both satellite and geophysical ERB echoes appear at the same fre- 
quency shift in the power spectra. While this feature renders the field-aligned pointing 
direction ineffective for a statistical discrimination approach, it actually provides un- 
ambiguous geophysical ERB identification when events are considered individually 
for the following reason. 

Enhanced echoes of geophysical origin were believed to be associated with intense 
electron drifts or currents [FOS88]. The most intense currents are also field-aligned 
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Figure 6-1: Histograms of velocity distribution of a) observed peaks, and b) satellite 
population. At this look angle geophysical enhancements and polar orbiting satellites 
share the same velocity range, [after Schlegel and Moorcroft, 1989]. 
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Figure 6-2: Histograms of altitude distribution of a) all observed coherent peaks, and 
b) satellite population. Note large number of observed peaks below 300 km. [after 
Schlegel and Moorcroft, 1989]. 
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and extend for many kilometers along the magnetic field direction. Enhanced echoes 
driven by these currents are therefore stimulated over the same extent along the 
field lines and can be observed with a radar pointed parallel to the same line. (The 
altitude extent of the Millstone Hill zenith observations is attributed to this effect; see 
Figure 5-17). Range extended echoes were discovered by RD91 when investigating 
the details of a cluster of spectra which caused correlator errors due to overflows 
in the A/D converters of the radar receiver. Two experiments were analyzed (Feb. 
14, 1990; Oct. 25, 1989), in which 24 ten-second integration periods were found to 
contain enhanced echoes. The following summary of those observations is adapted 
from RD91 and private communication with M. Rietveld, 1991. 

A sample of the data is shown in Figure 6-3, where five consecutive 10-sec in- 
tegration periods are shown which illustrate the transition from normal incoherent 
spectra to enhanced spectra and back again. The enhanced echoes are extended 
over a wide range of altitudes; a satellite penetrating the radar beam would appear 
in at most three range gates in the sampling technique used here. The frequency 
shift of the anomalous spectra also increases smoothly with altitude, corresponding 
to the increase in the ion acoustic speed with height, further evidence of the echoes 
geophysical origins. 

The spectra in Figure 6-3 vary dramatically during the 30 seconds in which the 
anomalous echoes appear. The spectra are enhanced at altitudes as low as 130 km 
and as high as 600 km by as much as 17 dB during the 04:49:50-04:50:00 period. 
The enhancements can be very asymmetric in either the upshifted or downshifted 
peak, or relatively symmetric with both ion acoustic peaks amplified by nearly the 
same amount. A transition in the sign of the frequency shift of the asymmetric peak 
enhancements can occur within a single 10 second integration period over several 
tens of kilometers. The statistics of which type of enhancement (asymmetric upshift, 
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Figure 6-3: Five consecutive 10-sec integration periods observed with the EISCAT 
UHF radar. Within the 50-sec time period the ion lines change from normal to 
enhanced and back again [after RD91]. 
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Figure 6-4: Histograms showing the height distribution of enhanced ion acoustic 
spectra categorized by a) total enhanced spectra; b) upshifted or d ownshifted peak 
enhancement ; and c) both peaks enhanced symmetrically [after RD91]. 


asymmetric downshift, or symmetric) was observed in the two experiments under 
discussion are presented in histograms in Figure 6-4. 

The figure indicates that there appear to be two main altitude regimes favoring the 
spectral enhancements, with a minimum rate of occurrence at around 250 km. In the 
higher altitude regime, the enhancements are primarily asymmetrically downshifted 
or symmetric; at lower altitudes upshifted asymmetric spectra are more common. The 
enhancements are found to be largest in the topside ionosphere above 250 km, and 
have lasted as long as 50 seconds continuously; shorter sporadic bursts distributed 
over a period of a few tens of minutes are seen more frequently. As more anomalous 
echoes from other experiments have been found, it appears that the low altitude 
enhancements are rare. Additionally, several ERB events have been observed when 
the radar was not pointed field-aligned, at was directed at aspect angles up to 20°. 
The echoes are then found to be localized rather than extended in range, as in the 
observations reported from Millstone Hill. 
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Geophysical Conditions 

Geophysically, the events occur in the presence of sharp fluctuations in the local 
magnetic field (~ 500 nT N-S), auroral particle precipitation evidenced by airglow, 
elevated electron temperatures ( T e up to 8000°K), and electron-to-ion temperature 
ratios greater than two, though not extremely high (typically, 2 < T e /Ti < 3). Several 
more ERB events have been discovered in EISCAT data since RD91; some of these 
cases were more closely correlated with elevated ion temperatures, strong perpen- 
dicular electric fields, and little, if any, auroral precipitation [Wahlund et al . , 1991]. 
Based on an evaluation of the sparse existing data, the physical mechanisms respon- 
sible for the observed enhancements are explored in Chapter 7. A recently conducted 
experiment to obtain better data is described in the next section. 

6.1.3 Dedicated ERB Experiment Performed at EISCAT 

In late November 1990, EISCAT radar time to acquire new data was provided by Dr. 
Kristian Schlegel of the Max Planck Institute for Aeronomy, and a single experiment 
designed to measure ERB with the highest possible resolution in time and space was 
conducted. The methods and results of the experiment are presented below. 

Experimental Description 

Experimental considerations for conclusively obtaining good quality enhanced echoes 
were determined by our understanding of the phenomenon. The experiment was 
designed to provide high resolution (in frequency and time) power spectra over a 
large altitude extent, primarily at field aligned and vertical look angles. The field- 
aligned position facilitates observing enhancements over a large altitude extent, while 
the vertical position constrains the possible line-of-sight satellite velocities to subsonic 
speeds in the ionosphere and shares common scattering volume with the VHF radar. 
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experimental Description 


SP-GE-SAT 

11/29/90, 16:00 UT - 11/ 30/90, 08:00 UT 


Radar Z&mnmn Lflk&k 1) 

- pulse length 500 

- sampling interval 10 ju 

- gate separation 75 km 

- number of scatter gates 13 

- total range 225-1175 km 

Scan Cycle 

- field-aligned 15 min 

- vertical 5 min 

Channels 5 and 6 were used to monitor plasma lines. All other channels were used to 
record the ion line. 


Radar Faramcvcn 1 Block 21 

- Same as CP-3 

Staa.Cy.fii.fi 

- 20 min elev. scan, 11 points (abbreviated version of CP-3 elev. scan) 

± 28° elev. scan in north- south plane through Tromse and Kiruna 

Block 2 run approximately every two hours for background diagnostics. The VHF 
antenna was pointed vertical throughout the experiment. 

Figure 6-5: Radar operating parameters employed in experiment to investigate ERB 
at EISCAT on November 30, 1990. 

Ideally we hoped to obtain simultaneous UHF and VHF measurements for radar cross 
section measurements at two frequencies. 

A 20-minute elevation scan was performed approximately every two hours for 
diagnosis of background ionospheric parameters. Two plasma line receiver channels 
were also recorded to deduce ionospheric currents directly. A complete description of 
the radar parameters and scan cycles employed in the experiment is given in Figure 

6-5. 

To distinguish the enhanced echoes from actual hard targets we relied on the 
multi-static capability of the EISCAT UHF system to obtain full 3-D Doppler infor- 
mation. Coherent echoes detected within the common scattering volume which did 
not exhibit velocities consistent with hard targets could be attributed to geophysical 
effects. In addition to hard target discrimination, the tri-static system can provide 
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the aspect angle dependence of the enhanced radar cross section, an important clue 
to understanding the physics of this phenomenon. To better facilitate the task of ac- 
quiring coherent echoes in the common scattering volume, a “smart” real-time radar 
control program was employed. 

This control program, named MONITOR, was written by Anthony van Eyken uti- 
lizing an “ERB detection” algorithm developed by Keith Groves and Mike Rietveld. 
The received ion line spectra were monitored in real-time to detect enhancements. 
When enhanced spectra were detected, the remote receiving antennas were auto- 
matically repositioned to observe the altitude of maximum enhancement along the 
transmitted radar beam, as depicted in Figure 6-6. This process was usually com- 
pleted within 10 seconds, considerably less than the time scales over which the echoes 
can persist (a few tens of seconds). The program was tested by looking at satellites 
to simulate ERB. An additional valuable feature of MONITOR is the capability to 
automatically reduce the integration period when enhanced power levels are detected, 
acquiring high time resolution measurements with adequate signal-to-noise. 

Initially, the 30-hour allotment of radar time was scheduled as a single time block 
beginning at 08:00 UT on November 29, 1990. Because the most promising periods 
for enhanced echo observations are believed to begin in the early evening and continue 
into the early morning hours, however, we chose to use only half the allotted time 
for the first experiment, beginning at 16:00 UT, and save the remaining time for 
additional observations on the evening of December 1. 

From a technical standpoint, the experiment was executed smoothly; unfortu- 
nately, it appears that enhanced radar backscatter was not observed. This is un- 
doubtedly due to the lack of appropriate geophysical conditions during the entire 
observational period. Magnetically, the ionosphere was very quiet; the magnetogram 
at the Tromsp site showed only weak, gradual variations, and the Ap index for the 
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Figure 6-6: Conceptual drawing of the ERB experiment performed at EISCAT in 
November, 1990. A realtime monitor program was employed to detect and locate 
ERB; the remote receiving antennas were then repositioned to sample the same vol- 


umi 


period was a meager five. Also, there was no evidence of particle precipitation, and 
electron temperatures remained below 3000° K. Additional observing periods were 
tentatively planned for the evenings of December 1 and 4. On the basis of the latest 
forecasts for very quiet conditions on both days, however, it was decided not to run 
dedicated coherent echoes experiments. A search through existing data was believed 
a more efficient means of locating additional ERB events. 


6.1.4 UHF Common Program Data 

The computer algorithm used to detect enhanced power spectra in real-time can also 
be used to analyze existing raw data tapes. Each spectra is checked, and those that 
exceed predetermined levels of asymmetry and/or overall power enhancement are 
recorded for later analysis. For these analyses, a 3 dB asymmetry threshold was used 
to reduce the large number of false detections due to low SNR, even though using 
the algorithm at this sensitivity level probably results in failure to detect some weak 
enhancements. I decided to look at CP-2 data initially, hoping to find enhanced UHF 
echoes in the vertical antenna position which could be compared with the VHF returns 
from the same scattering volume. CP-2 designates Common Program 2, a standard 
operation mode which includes both field-aligned and vertical radar look-angles. A 
total of 111 hours of Tromsq CP-2 data and 318 hours of data from the remote sites 
was reduced and analyzed for additional evidence of enhanced radar backscatter. The 
actual experiments analyzed are shown in Table 6.1. 

Despite the relatively large amount of data analyzed, significant occurrences of 
ERB were not found. As with the experiment, the probable cause for the lack of 
enhanced echoes is the low level of magnetic activity present during the time periods 
analyzed. In addition, only a small volume of space is monitored by the stationary 
radar beam which requires that ERB processes occur in very specific locations. While 
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Table 6.1: Table of Analyzed Common Program Experiments. T, K, and S denote 
Tromsp, Kir una, and Sodankyla. 


STRT. DATE-UT 

END DATE-UT 

MODE 

SITES 

87/10/20-09:30 

87/10/21-14:00 

CP-2-D 

K,S 

87/10/21-14:00 

87/10/22-23:30 

CP-2-D 

T,K,S 

88/08/09-08:00 

88/08/10-08:30 

CP-2-D 

K,S 

88/08/16-08:00 

88/08/18-23:00 

CP-2-D 

T,K,S 

89/10/23-18:00 

89/10/25-14:30 

CP-2-D 

K,S 

89/10/24-17:00 

89/10/25-07:30 

CP-2-D 

T 

90/11/27-17:00 

90/11/27-18:00 

CP-2-D 

T 


not indicative of all localized auroral processes, the Ap index for magnetic activity 
does provide an estimate of the general level of magnetic disturbance globally. This 
index was less than 15 (quiet) on all but one of the nine days investigated, and less 
them 10 on all but two. These periods were chosen for analysis initially because they 
were available at the Tromsp site. Later, due to several reasons (primarily computer 
failure), it was possible to process only half the data originally slated for analysis at 
EISCAT Headquarters in Kiruna, Sweden. In addition, 24 hours of data processed 
from disturbed time periods was lost due to a faulty data transfer. 


Tromso Data 

About 2.5% of the Tromsp spectra surveyed were determined to be asymmetrically 
enhanced; more than 95% of the selected spectra are not geophysically significant. 
Most result from satellites penetrating the radar beam; some are asymmetric because 
of random noise spikes in low SNR data from high altitudes. Only two cases of interest 
were found in the data examined thus far. 

The first, shown in Figure 6-7, is comprised of enhanced spectra previously found 
by Mike Rietveld and proves the value of this method of analysis. The data quality for 
the October 24/25 experiment was quite good, and only 42 spectra were computer- 
selected from 14 hours of data. Of these 42, 16 are clearly enhanced by geophysical 
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Figure 6-7: Asymmetric enhanced UHF spectra from October 1990 extended over 13 
consecutive range gates. The range gate number is on the upper left of each plot; 
the number on the upper right is the asymmetry index assigned by the MONITOR 
detection algorithm. 

mechanisms as evidenced by their sequential extent in range. This data set is the 
exception with regard to geomagnetic activity noted above. It was recorded three 
days after the peak of one of the largest magnetic storms on record. Actual conditions 
for the day and hour of the enhanced echoes occurrence are characterized by an Ap 
index of 23 and a 3-hr Kp of 5. 

Data from August 18, 1988 form the second interesting case (see Figure 6-8), when 
weakly enhanced spectra were observed in the UHF radar beam continuously for more 
than 30 seconds at ranges varying from 520 - 600 km along the magnetic field. The 
localized nature of the enhancements does allow the possibility that the spectra may 
be a hard target signature in the beam. Figure 6-9, however, indicates that a satellite 
at 600 km range could remain in the stationary radar beam for no more than a few 
seconds, and direct measurements of satellites with the EISCAT system support this 
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Figure 6-8: Weakly enhanced ion acoustic spectra seen in four consecutive 10-second 
integration periods. The range varies from 520 km to 600 km; a satellite at this range 
would pass through the radar beam in less than five seconds. 


conclusion. 

These spectra were selected by the MONITOR program and represent the inter- 
esting data I have from this time period. It is entirely possible that there are other 
spectra recorded during this time period which, while not detected by MONITOR 
(and therefore not recorded), contain meaningful levels of asymmetric enhancement. 
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Hard Target Resident Time in 
EISCAT UHF Radar Main Beam 



Figure 6-9: Approximate main beam satellite transit time as a function of altitude 
for the 1° EISCAT UHF radar beam. Targets below 1000 km remain in the beam for 
less than five seconds. Circular orbits are assumed. 

Remote Sites Data 

As pointed out previously, enhanced echoes observed by all sites can be conclusively 
attributed to geophysical sources if they do not exhibit a 3-dimensional Doppler ve- 
locity consistent with orbiting bodies. Motivated by the promise of unambiguous 
classification, a large amount of data from the remote receiving sites at Kiruna and 
Sodankyla was examined. However, not a single case of ERB was detected by MON- 
ITOR in more than 300 hours of remote site data. 

The lack of detections is somewhat puzzling, particularly because during the Oc- 
tober 24/25 experiment moderately enhanced spectra were observed at Tromsp from 
the common volume scattering height of 278 km. Unfortunately, the remote sites were 
not pointed at this altitude during the period of observed ERB; in fact, no coincident 
measurements of ERB have been discovered thus far. This can be partially explained 
by the fact that the remote sites monitor altitudes at 278 km or below, where very 
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few enhancements have been seen. 

The aspect angle sensitivity of ERB processes may prove important, as well, be- 
cause the remote sites are not able to view the F region common scattering volume 
at small aspect angles. Because the enhancement phenomenon is associated with 
field-aligned processes, it is possible that it produces a much weaker signature when 
diagnosed by the remote sites at large aspect angles. A final consideration is that 
the scattering volume available to the remote sites is determined by the intersection 
of their receiving antenna beam and the radar beam transmitted from Troms0. The 
respective radax beams are only about five kilometers in diameter at the intersec- 
tion altitude of 275 km. This provides a relatively small scattering volume, and the 
power received by the remote sites is considerably less than that received from direct 
backscatter at Troms0. Nevertheless, the continuing effort to find more cases of ERB 
in the EISCAT data will probably uncover some events diagnosed simultaneously by 
both the Troms0 receiver and the remote sites, providing new insights into the physics 
of these events as discussed in Chapter 7 . 

6.1.5 EISCAT VHF Observations of ERB 

Evidence of ERB events detected with the EISCAT VHF radar were recently reported 
by Collis et al., [1991] (henceforth, C0LL91). The enhanced spectra are similar in 
appearance to those recorded by the UHF radar discussed previously. A summary of 
their characteristics follows. 

Summary of VHF Observations 

VHF ERB was observed during experiments in January 1989 and February 1990 
[COLL91]. Both cases were associated with unusually strong visible red aurora as 
indicated by all-sky imagers and scanning photometers. The auroral signatures are 
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believed to be caused by low energy ( E < 100 eV) precipitating electrons. A sam- 
ple of the VHF echoes is shown in Figure 6-10, where the altitude variation of the 
enhancements over a three minute period in the February experiment is evident. 

The apparent downward motion of the ERB is interpreted by COLL91 to represent 
the equatorward motion of a field-aligned structure ( e.g ., auroral arc) passing over the 
radar. The VHF echoes shown in Figure 6-10 are not extended in range because the 
antenna is pointed to the zenith, rather than field-aligned as in the case of the UHF 
enhancements. Further consideration of the data indicates other differences between 
the VHF and UHF enhancements. The VHF events are found at much higher altitudes 
(up to 1350 km) and persist for much longer, up to four minutes in the example shown 
here. In the February 1990 experiment VLF spectra were enhanced sporadically 
for more than one hour. By comparison, UHF events at EISCAT have not been 
found above 600 km, and 50 seconds represents the longest period of uninterrupted 
enhancement, while the duration of sporadic enhancement periods does not exceed 
about ten minutes. These facts indicate that the enhancements may be more easily 
excited at VHF frequencies [COLL91]. 

The frequency shift of the enhanced peak was seen to change from upshifted to 
downshifted during the course of about a minute in the January 1989 experiment. 
ERB data acquired during the February 1990 event, however, are characterized by 
a consistent dominant enhancement of the downshifted peak for up to four minutes, 
though at times both peaks were enhanced in the same integration period. This is 
in contrast to the UHF measurements where the asymmetry seems to fluctuate more 
rapidly in both time and space between downshifted and upshifted frequencies. 

The total cross section enhancement for the VHF spectra ranges from factors 
of 4-5 during the January experiment to greater than 25 in the February events. 
Measurements of ERB by both the UHF and VHF radars have been recorded during 


129 




-4 0 *4 -4 0 +4 -4 0*4 

Frequency (kHz) 


Figure 6-10: Enhanced zenith VHF spectra from 20 February, 1990. The three 10- 
second integrations shown were recorded 1 minute apart beginning at 18:03:00 UT; the 
altitude of the enhancement varies smoothly. Enhancements above normal incoherent 
scatter are shaded black [after COLL91]. 

the same time period, confirming that both are associated with elevated electron 
temperatures, large magnetic field fluctuations, and, when optical data is available, 
with red aurora. Unfortunately, comparisons of simultaneous UHF/VHF spectra 
from the same scattering volume have not been reported. This is because the VHF 
antenna is nearly always pointed vertical during experiments, while the UHF dish 
acquires data at other look-angles a majority of the observing time. The comparison 
of spectral characteristics at two frequencies would provide important information on 
the scale size dependence of the physical mechanism responsible for the ERB. 

8.2 ERB Observations at Sondrestrom 

The incoherent radar located at Sondrestromfjord, Greenland operates at 1290 MHz 
and offers an ideal location («74° Inv.lat.) for the excitation and detection of ERB at 
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a wavelength shorter than that of both the EISCAT and Millstone Hill systems. The 
data base available for such studies is, at present, much more limited than either of 
the latter sites. Frequent low energy electron precipitation occurring at this latitude 
increases the likelihood of backscatter enhancements locally. 


6.2.1 Zenith Experiments at Sondrestom 

As explained previously, acquiring data with the radar antenna pointed to zenith 
effectively filters out hard target contamination by limiting the line-of-sight velocity to 
subsonic magnitudes. Accordingly, attempts to identify ERB events at Sondrestrom 
have employed this approach. 

Analysis of four hours of zenith data taken in August, 1988 revealed that two of 
the 20 coherent echoes observed occurred at frequency shifts corresponding to ion 
acoustic waves; the remainder of the enhanced peaks had essentially no Doppler shift 
and are presumed to be caused by hard targets passing through the radar beam. The 
two events at the ion acoustic frequency showed backscattered power enhanced by 
factors of 2 and 5, respectively; smaller asymmetries were evident in several other 
spectra. The spectrum enhanced by a factor of 5 is shown in Figure 6-11. The range 
gates shown in the figure do not overlap, so that a localized target appears in only 
one gate. This is consistent with EISCAT and Millstone Hill observations of ERB for 
nonzero aspect angles. 

Six more zenith experiments were run at Sondrestrom in May, September, and Oc- 
tober, 1989 and October, 1990 providing 44 additional hours of observations. These 
data, however, show little evidence of ERB at ion acoustic frequencies; asymmetric 
enhancements are at most 1.5 times greater than the normal background spectra, 
not considered statistically significant. No extended periods of ERB occurrence have 
been found in Sondrestrom data yet. The lack of ion acoustic ERB in these ex- 
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Figure 6-11: Sondrestrom zenith data showing factor of 5 enhancement in ion acoustic 
spectrum. Altitude scale is in kilometers. 




periments suggests the absence of intense field-aligned currents. The magnitude of 
these currents is related to the level of geomagnetic activity, generally weak during 
these experiments. Another reason for the absence of geophysical ERB may be that 
higher electron drifts are required for enhancement at the short wavelength (12 cm) 
diagnosed by the Sondrestrom radar. 


6.3 Summary of ERB from Other Sites 

Important progress has been achieved towards identifying and understanding en- 
hanced radar backscatter from the ionosphere. Such events are known to occur briefly 
and infrequently under special ionospheric conditions; enhancements below 250 km 
are especially rare. The most prominent observations reported by RD91 and COLL91 
associate ERB events with red aurora (viz., particle precipitation), large magnetic 
field variations locally, and elevated electron temperatures. These requirements are 
also reflected by global magnetic indices for the periods indicating very disturbed 
conditions. The EISCAT ERB observations have been seen when Kp> 5, a level of 
disturbance occurring about five percent of the time during solar max [ Ma.yo.ud, 1980]. 
Comparing the Millstone Hill, EISCAT and Sondrestrom observations suggests that 
ion acoustic spectra may be enhanced over a wide range of scale sizes; enhancement 
thresholds appear to decrease with increasing wavelength. 

These findings are basically consistent with the current driven source mechanisms 
of Kindel and Kennel [1971] and Rosenbluth and Rostoker [1962] suggested by FOS88, 
RD91, and COLL91. These mechanisms show that a relative drift between the elec- 
trons and the ions ( i.e ., a current) in the plasma will produce asymmetric enhance- 
ments in the observed ion acoustic spectrum. The frequency shift of the asymmetry 
depends on whether the electrons are drifting towards (upshift) the radar or away 
(downshift) from the radar. A puzzling inconsistency with the known mechanisms, 
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however, stems from the commonly observed simultaneous enhancement of both up- 
shifted and downshifted ion acoustic peaks. The magnitude of enhancements larger 
than 7 dB and the detection of ERB at nonzero aspect angles axe also not described 
adequately by RR and KK. Possible explanations for these observations based on 
existing theories and an alternative mechanism for ion acoustic enhancements are 
presented in Chapter 7. 
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